DUDLEY  Kvrn-'  LIBRARY 

)UATE  SCHOOL 
..  CALIFORNIA  93943 


NAVAL  POSTGRADUATE  SC 

Monterey,  California 


THESIS 

COMBAT 

ATTRITION  ANALYSIS  USING  RENEWAI 
by 

Cho ,  Joong  Kun 
March  1984 

i  PROCESS 

Thesis  Advisor:                      J.  G. 

Taylor 

Approved  for  public  release;  distribution  unlimited 


UNCLASSIFIED 


SECURITY  CLASSIFICATION   OF   THIS  PACE   (Wh»n  Data   Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.    REPORT  NUMBER 


2.  GOVT  ACCESSION   NO 


3.     RECIPIENT'S  CATALOG   NUMBER 


4.    TITLE  (and  Submit) 


Combat  Attrition  Analysis  Using  Renewal  Process 


5.     TYPE  OF   REPORT  a   PERIOO  COVERED 

Master's  Thesis 
March,  1984 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS 

Cho,  Joong  Kun 


8.  CONTRACT  OR  GRANT  NUMBERf»J 


9.  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


10.     PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4    WORK  UNIT   NUMBERS 


II.     CONTROLLING  OFFICE  NAME  AND  AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93943 

14.     MONITORING  AGENCY  NAME   4    AOORESSf//  different  from  Controlling  Office) 


12.     REPORT   DATE 

March,    1984 


13.     NUMBER  OF  PAGES 
55 


15.     SECURITY   CLASS,   (at  this  report) 

Unclassified 


15«.     DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  (of  thla  Report; 


Approved  for  Public  release;  distrubution  unlimited 


17.     DISTRIBUTION  STATEMENT  (of  the  eburaet  entarad  In  Block  20,   It  different  from  Report) 


18.     SUPPLEMENTARY  NOTES 


19-     KEY  WORDS     Continue  on  reveree  aid*  It  naeaaaary  and  Identity  by  block  number) 

Serial  Acquisition  Renewal  Density 

Parallel  Acquisition 
Renewal  Function 
Renewal  Ecuation 


20.     ABSTRACT  (Contlnua  an  ravaraa  tide  It  naeaaaary  and  tdantlty  by  block  number) 

This  thesis  uses  renewal  theory  to  investigate  the  Lanchester-type 
combat  attrition  process.    The  attrition  process  is  analysed  in  detail 
and  modelled  as  a  so-called  renewal  process  in  which  times  between 
casualties  are  considered  to  be  independently  and  identically  distri- 
buted random  variables.    Other  random  variables  that  can  be  considered 
in  the  renewal  process  are  examined,  and  the  distrubutions  of  these 
random  variables  are  determined  in  order  to  study  the  behaviour  of  attri- 


DD  ,  '™n  1473 


EDITION  OF   1  NOV  SS  IS  OBSOLETE 
S-'N  0102-  LF-OU-6601 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF   THIS  PAGE  (Whan  Data  Sntare. 


yiNvgij£oaxj:  lh,u 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whin  Data  Entarad) 


#20.  ABSTRACT  (Continued) 

tion  process.    Examples  with  specific  distribution  functions  are 
given  for  better  understanding.    Computer  simulation  is  generated 
and  compared  with  the  attrition  model  developed.    The  total  casualty 
occurence  by  total  force  is  also  discussed,  through  pooling  of  the 
single  renewal  processes.    The  total  casualty  occurence  is  shown  to 
be  a  poisson  process  and  times  between  casualties  to  be  approximately 
exponentially  distributed  for  large  numbers  of  combatants. 


S  N  0102-  LF-  014-  6601 


UNCLASSIFIED 

2       SECURITY   CLASSIFICATION  OF  THIS  P  AGE(V*>»n  Data  Entarad) 


Approved  for   public  release;    distribution    unlimited. 

Combat  Attrition  Analysis   Using    Renewal   Process 

by 

Joon g   Kun,Cho 

Major,   Reoub lie   of   Korea    Army 

B.A.,    Korean   Military   Academy , 1 974 


Submitted    in   partial   fulfillment   of    the 
requirements    for   the   degree    of 


MASTER    OF    SCIENCE    IN    OPERATIONS    RESEARCH 


from   the 


NAVAL  POSTGRADUATE  SCHOO! 
March  19  34 


Dl 

Nay: 

M0K  CA   9394; 


ABSTBACT 


This  thesis  uses  renewal  theory  to  investigate  the 
lanchester-type  comfcat  attrition  process.  The  attrition 
process  is  analysed  in  detail  and  modelled  as  a  so-called 
renewal  process  in  which  times  between  casualties  are 
considered  to  be  independently  and  identically  distributed 
random      variables.  Other      random    variables      that      can      be 

considered  in  the  renewal  process  are  examined,  and  the 
distributions  of  these  random  variables  are  determined  in 
order  to  study  the  behaviour  of  attrition  process.  Examples 
with  specific  distribution  functions  are  given  for  better 
understanding.  Computer  simulation  is  generated  and  compared 
with      the    attrition      model    developed.  The    total      casualty 

occurence  by  total  force  is  also  discussed,  through  pooling 
of  the  single  renewal  processes.  The  total  casualty  occu- 
rence is  shown  to  be  a  roisson  process  and  times  between 
casualties  to  be  approximately  exponentially  distributed  for 
large    numbers    of   combatants. 
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I.    INTRODUCTION 

Today/  deter ministic-dif f erential-equation  models  (  that 
are  commonly  called  Lanchester-type  combat  models)  are 
widely  used  in  military,  operations  research  for  defense 
planning  purposes,  especially  as  concerns  land  combat.  In 
combat  analysis,  such  questions  frequently  arise  as  'Who 
will  win  the  battle  ?,  How  long  will  the  battle  last  ?,  Kow 
many  survivors  will  the  winner  have  ?,  and  Kow  do  the  force 
levels  change  over  time  ?',etc.  There  have  neen  a  lot  of 
studies  to  estimate  attrition  in  battle  to  answer  such  ques- 
tions  as   these. 

In  1914 ,F.  W. Lanchester  [Ref.  1]  first  proposed  a  mathe- 
matical model  of  combat  in  order  to  justify  the  principle  of 
concentration  of  forces  under  "modern"  conditions.  He 
believed  that  in  "ancient"  times  war  was  a  series  of 
one-to-one  duels  between  men  engaged  in  hand-to-hand  combat 
using  sword,  axe,  lance,  etc.  So  there  would  be  no  advan- 
tage gained  from  concentration  of  forces.  Lanchester 
reasoned  that  modern  technology  had  changed  the  nature  of 
this  "ancient"  warfare  and  now  there  was  a  decided  advantage 
to  "concentration  of  forces".1  He  hypothized  that  each 
side's  casualty  rate  was  proportional  to  the  number  of 
opposing  firers.  Implicit  (  but  not  yet  explicit)  in 
Lanchester1 s  original  work  is  the  concept  of  an  "attrition- 
rate  coefficient"  as  the  rate  at  which  a  single  firer  kills 
a  particular     enemy    target    type.           Furthermore,      Lanchester 


1  Helmtold  notes  that  the  analysis  of  the  data  from  92 
land  combat  battles  suggests  that  victory  in  cattle  is 
primarily  determined  by  the  factors  other  than  numerical 
superiority  and  challenges  the  ability  of  any  model  of 
combat  which  concentrates  almost  exclusively  on  numerical 
force  size  to  yield  a  practically  useful  predictor  of 
victory    in    battle. 


thought  that  the  number  of  survivors  were  the  critical 
factor  deciding  the  outcome  of  the  battle  in  "modern" 
warfare. 

Since  Lanchester's  original  efforts  to  describe  the 
dynamics  of  combat  mathematically,  many  analysts  -  in  an 
attempt  to  add  reality  to  the  combat  description  -  have 
extended  the  theory  to  include  additional  factors.  Many  of 
these  extentions  were  described  by  Dolansky  [Ref.  2].  In 
his  paper,  Dolansky  pointed  out  that  the  Lanchester 
attrition-rate  coefficients  were  nard  to  determine  for 
particular  weapons,  and,  accordingly,  coefficients  of  known 
magnitude  had  been  assumed  in  most  models.  Thus,  Dolansky 
concluded  that  the  use  of  Lanchester-type  models  for 
prediction  of  battle  results  had  been  hampered  by  this 
inability  to  predict  numerical  values  for  attrition-rate 
coefficients.  However,  even  if  we  may  be  in  doubt  as  to  the 
proper  form  of  attrition  relations,  we  may  at  least  be  able 
to  make  pragmatic  assertions  such  as,  'If  the  Lanchester 
attrition  is  used  for  predictions,  then  by  measuring  certain 
factors  and  perfQrming  certain  operations  a  prediction  of 
such-and-such  an  aspect  of  combat  can  be  obtained. '  So, 
although  there  has  been  a  continuing  discussion  among  mili- 
tary operations  analysts  about  the  merits  of  Lanchester-type 
attrition-rate,  it  is  still  generally  accepted  as  the  basic 
founding  for  all  the  combat  models  studies  and  attrition- 
rate    coefficient   is    still   a    core    of    these    studies. 

Bonder/Barf oot  [ Ref .  3]  proposed  that  an  attrition-rate 
coefficient  be  defined  as  the  reciprocal  of  the  expected 
time  for  an  individual  firer  to  kill  an  enemy  target.  A 
more  thorough  dicussion,  however,  of  the  justification  for 
Bonder/Barefoot  is  to  be  found  in  Taylor  [Ref.  4].  Barf  oot 
[Ref.  5]  suggested  that  the  attrition  process  be  considered 
as  a  renewal  process .  Thus,  for  a  single  Y  firer,  we  will 
assume   that   each      individual   Y   kills    X      targets    according    to 
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an  attrition  process  in  which  times  to  kills  are  indepen- 
dently and  identically  distributed  (i.i.d.)  random  variables 
(r.v).  Vie  will  also  assume  that  each  firer  acts  indepen- 
dently of  any  other  firer,  and  similarly  for  targets.  Such 
an  attrition  process  turns  out  to  be  a  renewal  process  in 
the    theory   of   stochastic  processes. 

Once   we  introduce   the    theory   of   renewal   process,    we   can 
consider   many    other    random    variables    such    as    the    time    up    to 

n-th    kill,    S    ;       the      number    of     .kills      in      time    t,    N  (t)  ;    the 

n 
expected   number    of    kills   in    time      t,      :1  (t)  ,      etc.      Thus   this 

thesis  starts  with  consideration  of  3onder/Taylor 's  concept 
of  attrition-cat e  coefficient,  and  then  studies  the  prob- 
ability distribution  of  the  time  between  casualties 
according   to  the  theory   of    renewal   process. 

In  this  thesis,  we  will  mainly  focus  on  the  fundamental 
Lanchester-type  attrition  paradigm  which  deals  with  homoge- 
neous forces,  because  understanding  this  basic  paradigm  is 
essential  for   extentions   to    further    complex   models. 

This  thesis  is  organized  in  the  following  fashion. 
First,  we  review  a  few  relevant  facts  about  the  Lanchester 
attrition  models.  Then  we  review  basic  ideas  of  renewal 
process  and  develop  the  random  variables  whose  study  is  of 
our  concern,  followed  by  the  study  of  these  random  vari- 
ables. Then  we  will  describe  the  combat  attrition  process  by 
comparing  the  values  of  different  random  variables  obtained 
through  renewal  process.  Finally,  the  conclusion  mainly 
contributes  to  the  summary  of  the  work  and  suggests  further 
studies. 
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II.    FUNDAMENTAL    LANCHESTER-TYPE    ATTRITION    PARADIGM 

A.       THE    EASIC    PARADIGM 

Let  us  consider  combat  between  two  homogeneous  forces:  a 
homogeneous  X  force  opposed  by  a  homogeneous  Y  force.  We 
will  focus  on  the  force-on -force  attrition  process  in  the 
combat  between  these  two  homogeneous  forces.  The  basic 
Lanchester-type  paradigm  for  "modern"  warfare  assumes  that 
the  casualty  rate  of  such  a  homogeneous  force  is  directly 
proportional  to  the  number  of  enemy  forces .  e.g.  the  X  force 
casualty   rate    is   given  by 

dx 

—    =    -    ay  with        x(0)     =    x      , 

dt  0 

and    the   Y   for^e    casualty  rate   is    given   by 

ij 

-=■    =    -    bx  with        y(0)     =    y      , 

dt  0 

where    ta'    denotes    the  rate    at    which    a    single    Y   firer   kills    X 
targets   and   is   called  a   Lan cnester   at tr it ion- rate    coef f i- 
cient.    Here   x      and    y  denote        the         number      of      X      and      Y 

comcatants  (respectively)  at  the  beginning  of  tne  battle. 
Many  factors  may  affect  the  value  of  'a1  and  '  b',  but  for 
the  time  being  it  is  not  essential  that  we  be  explicit  about 
functional  dependence  of  a  and  b.  Let  us  next  address  the 
problem  of  computing  a  reliable  numerical  values  for  such 
attrition-rate    ccef f icients . 

Two  approaches  have  been  developed  in  the  United  States 
to    determine  the   value   of    these   coefficients    as    follows: 

(1)         A  statistical      estimate      based      on    "combat"      data 
generated    by   a    detailed    Monte   Carlo    combat    simulation. 
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(2)  An  analytical  submodel  of  attrition  process  for  the 
particular  combination  of  firer  and  target  type. 
In  this  thesis  we  will  consider  only  the  second  approach. 
S. Bonder  [Ref.  6]  has  called  this  approach  the  use  of  a 
free-standing  or  independent  analytical  model.  The  basic 
idea  of  this  approach  is  to  develop  an  analytical  expression 
for  each  required  kill  rate  by  considering  the  single  firer 
engaging  a  'passive'  target  (i.e.  one  that  doesn't  fire 
back)  and  then  tie  all  the  attrition  rates  together  in 
f orce-on-f orce    combat   with    a   Lanchester-type    model. 

Barfoot  [Ref.  5]  in  his  study  has  suggested  taking  a 
Lan Chester  attrition-rate  coefficient  as  the  reciprocal  of 
the  ex£ected  time  for  an  i  n dividual  firer  to  kill  an  en  em  y 
'target.  '2  TJithin  context  of  the  above  homogeneous  force 
Lanchester-type    combat   model, this    means 

1 

E[*T       ]    , 
XY 

where    T        denotes   a    r.v.    representing    the    time   for    a    single 

XY 
Y   firer    to      kill   an    X   target   and    E[T]      denotes   expected   time 

to    kill.      Taylor   also  supported    Barfoot's   idea   in    his   recent 

study   [Ref.    7]      by    comparing   times      between    Kills      that   were 

both        exponentially        distributed        and        non-exponentially 

distributed,      and      using   previous    equation    to      determine   the 

expected   time   to  Kill  a    target. 


2However,  the  justification  of  this  approach  is  not 
accepted  and  is  apparently  somewhat  controversial.  But  in 
case  of  exponentially  distributed  times  between  kills,  a 
deterministic  Lanchester-type  model  may  indeed  be  considered 
to  yield  the  mean  course  of  comnat.  Anyway,  this  approach  is 
accepted  as  general  guide  for  further  study  of  attrition- 
rate    coefficient. 
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B.       DETEBMINATION   OF    EXPECTED    TIME    TO    KILL 

Bonder  and  Farrel  £Eef.  8]  have  developed  general  meth- 
odology for  determining  the  attrition-rate  coefficient  by 
determining  the  expected  time  to  kill  a  target,  E[T],  for  a 
wide  spectrum  of  weapon  system  types.  Morever,  research 
since  the  mid-1960's  has  led  to  the  development  of  several 
other  methods  for  computing  the  expected  time  to  kill  a 
target.  For  present  purposes,  two  methods  will  be  of  our 
interest, 

(1)  Method    based  on   sum    of    component    event    times, 
and 

(2)  Method   based  on      first-passage   time   in    Semi-Markov 
process 

If  we  assume  the  special  case  of  tactical  interest:  namely, 
the  case  of  Mar koy-dependen t  fire, Bonder  [Ref-  3]  has  shown 
that    E[  T  ]  turns    out    to  ne  as   follows: 

(t   *t  )      (t    +t    ) 

h      r  m      i       i  1-P  (H|  H) 

Z£T]=t    +   t    -t    + * { +    P(H|H)     -    ?    } 

a         1       n       P(KJH)        P(H|M)  P(K|H)  1 

where  all  symbols  are  defined  in  Table  I  This  expession  for 
E[  T  ]    holds    for    the    following   assumptions; 

(1)  Markov-dependent    fire    with      parameter    ?    ,       P  (H|  H)  , 
and   P  (H|M) 

(2)  Geometric    distrioution      for      tne      number    of      hits 
required   for  a    kill   with    parameter   P(KjH) 

In  assessing  the  time  factor  to  determine  the 
attrition-rate  coefficient,  it  turned  out  that  target  acqui- 
sition process  and  weapon  system  capability  were  most  impor- 
tant factor  from  the  formula.  However,  wnen  we  consider  the 
time  factor  of  the  attrition  process,  target  acquisition 
process  would  be  the  most  important  factor  to  study  because 
it  directly  affects  the  required  length  of  time  to  kill  a 
target.  So  we  will  elaborate  on  target  acquisition  process 
for    further    use    in    attrition   analysis. 
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TABLE    I 
Variables  for  Expected   Time    to   Kill 

time    to   acquire   a    target,    t 

a 

time   to   fire    first  round    after   target   acquired,    t 

time    to    fire    a   round    following    a    hit,t 

h 

time   to   fire    a   round   following   a   miss,t 

m 

time   of   flight   of    the  projectile,    t 

probability    of  a    hit    on    first    round,    P 

1 

probability    of   a    hit    on    a   round   following    a   hit,    P(H|H) 
probability    of  a    Lit    on    a   round    following    a   miss,    ?  (H| M) 
probability    of   destroying    a    target    given   hit,    P  (K j H) 


C-       TARGET    ACQUISITION    PROCESS 

Target  acquisition  is  the  initial  step  in  engaging  a 
target  and  a  most  deciding  factor  that  affects  the  time  to 
kill.  An  important  distinction  made  in  VRI's  Lanchester-type 
combat  models  is  whether  the  target  acquisition  process  of  a 
single  "typical"  firer  type  is  a  serial  process  or  a 
parallel  p r oce s s  [Ref.  4].  The  two  modes  for  the  target 
acquisition  process  considered  by  VRI's  models (  including 
V2CT0R-2   )     are    as   follows: 

(1)  Serial   Acquisition 

(2)  Parallel    Acquisition 
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Here  a  single  firer  using  serial  acquisition  does  not 
acquire  targets  while  engaging  another  target,  say  target  A. 
When  such  a  firer  ceases  to  engage  a  target  A  (  due  to  kill, 
lost  or  any  other  reason),  then  he  must  acquire  a  new  target 
all  the  way  from  the  beginning-  So  it  is  assumed  that  he 
does  not  remember  any  acquisitions  made  prior  to  engaging 
the  target  A.  On  the  other  hand,  a  firer  using  parallel 
acquisition  searchs  for  targets  continuously,  'even  while 
engaging  target  A  and  remembers  tnose  targets  that  have  been 
acquired.  Thus  when  such  a  firer  finishes  the  engagement 
with  target  A,  he  can  immediately  shift  engagement  to  the 
next  target,  provided  chat  such  a  target  was  acquired  and 
processed  during  or  before  the  engagement  of  target  A.  For 
both  methods,  the  VRI  models  assume  that  a  firer  never 
engages    a   killed   target    again. 

1 •      Serial    Acquisition 

Now  when  we  consider  the  basic  Lanchester-type 
combat  model  in  the  previous  section,  a  firer  in  serial 
process  mode  of  target  acquisition  starts  again  from  the 
beginning  in  searching  a  new  target  whenever  the  previous 
target  has  teen  killed-  Thus  we  assume  that  target  acquisi- 
tion time  is  always  needed  to  acquire  a  different  target  for 
engagement-  Here  it  is  assumed  that  the  total-force  kill 
rate  is  just  the  sin gle-f ir er-kill-ra te  times  the  number  of 
firers,    e.g. 

dx 

"  it  =  ay  ' 

since  each  firer  is  independent  from  any  other  firer.  Also, 
in  the  case  of  a  serial  process  for  each  firer,  the  expected 
time  to  kill  a  target  is  the  same  as  has  been  shown  in  the 
equation  for  E[  T  ]  in  section  B.  A  summary  of  attrition-rate 
coefficient  results  for  serial  acquisition  is  given  in 
Table    II   . 
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TABLE   II 
Attrition   Process   Osing    Serial   Acquisition    Mode 


dx 

—   =   -   ay 

dt 

1 

E£T~~] 

XY 

T        =  time   for  a    Y   firer   to   kill    an    X    target 
XY 

(t    +t    )      (t    +t    ) 
h      f  id      f      £  l-p(HJH)  1 

E£I]    =    t   *t    -t    + ♦ £ +  ?(H|H)-P    } 

a       1       h    P(KJH)        P  (H|  M)  P  (Kj  H)  1 


2-      Carallel   Ac  ^  uisitio  n 

As  we  mentioned  earlier,  a  firer  in  the  parallel 
mode  of  acquisition  continues  to  acquire  new  targets  even 
while  he  is  engaging  a  particular  target.  Cnce  this  target 
has  been  killed,  then  he  car  immediately  shift  fire  to  a  new 
target      which   has      been    acquired      already.  In    this      case, 

total-force  kill  rate  is  given  by  the  product  of  the  kill 
rate  of  single  firer  against  acquired  targets  and  the 
expected  number  of  firers  who  have  already  acquired  one  or 
more    targets. 

The      summary      results   for      attrition      process      using 

parallel   acquisition   is   given   in    table    III    .    Note    here    that 

f        is   same   as    the    orcbabilitv    that   a    single      Y      firer    using 
xy 
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parallel   acquisition   mode   has    available    one   or  more    acquired 

targets    at    which   to    fire   at    time    t,       because    if    he    is    firing 

we  know    that  he    already   acquired    the    targets.    Then   it    should 

also   be   noted  that    f        v   represents   the    expected    number   of   Y 

xy 

firers  who  have  already  acquired  one  or  more  X  targets. 


TABLE  III 
Attrition  Process  Using  Parallel  Acquisition  Mode 


dx 

--  =  -  f  a'y 

dt      xy 


f      :    Prob(single    Y   firer    using    parallel   acquisition 
xy  is  firing   at    X    target    at    random  time    t        } 

<t 
f   =  1  -  exp  {  -  x    B   (s)  ds} 
xy  0   x  y 

1 

a*  =  :Y  firer  kill  rate  against  acquired  target 

E[T»       ] 
X  Y 


T'      =   Time    for    a    Y   firer   to   kill   an    acquired    X    target 
XY  (conditional   kill    time) 


E         :    the   rate    at   which    a    Y   firer    acquires   X    targets 
xy        at  time   t    when   appropriate 

(t  +t  )         (t  nj 

h      f  m      f       [  1  -  ?  ( H  J  H )  1 

E[T']    =    t    -t    + + ( +  P(H|H)-P    } 

1        h         ?(K|H)  P  (H|  3)  ?  (KjH)  1 
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Now  we  may  assume  that  it  would  take  longer  to  kill 
a  target  for  a  firer  using  serial  acquisition  than  a  firer 
using  parallel  acquisition  because  of  the  time  fraction  of 
acquisition  of  target.  This  distinction  between  these  two 
acquisition  modes  may  be  further  studied  as  the  case  of  a 
delayed  renewal  process.  Here  we  assume  that  expected  time 
to  kill  the  first  target  would  be  obtained  from  the  serial 
acquisition  but  following  time  to  kill  would  be  from  the 
parallel  acquisition,  provided  that  the  targets  are  acquired 
already.  But  in  this  thesis,  we  will  only  deal  with  the 
firer  who  uses  the  same  acquisition  mode  to  attrit  all 
targets. 

Before  continuing,  it  is  important  that  the  reader 
understand  both  the  basic  Lanchester-type  attrition  model 
and  also  Eonder's  approach  of  tne  attrition-rate  coefficient 
as  the  reciprocal  of  expected  time  to  Kill  a  target.  It  must 
be  recognized  that  this  expected  time  to  kill  a  target  is 
primarily   dependent    en   the    acquisition   process. 
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Ill-    ANALYSIS    OF    THE    COMBAT    ATTRITION    PHOCESS 

A-       INTRODUCTION 

The  Lanchester  attrition-rate  coefficient  is  the  rate  at 
which  a  single  firer  kills  a  particular  enemy  target  type  in 
lanchester-type    combat.  Development    of      technically-sound 

and  scientifically-valid  methodology  for  determining  numer- 
ical values  for  Lanchester  attrition-rate  coefficient  is  an 
essential  prerequisite  for  building  a  Lanchester-type  combat 
model.  The  basic  construct  cf  3onder/3arfoot  methodology  was 
to  take  a  Lanchester  attrition-rate  coefficient  as  the 
reciprocal  cf  the  expected  time  for  an  individual  firer  to 
kill  an  enemy  target  as  mentioned  in  the  previous  chapter. 
Taylor  provided  justification  [Ref.  4]  for  taking  Lanchester 
attrition-rate  coefficient  as  the  reciprocal  of  the  expected 
time  to  kill.  He  started  frcm  the  oasic  hypothesis  that 
combat  is  a  complex  random  process,  but  it  contains  enough 
regularity  that  the  appropriate  Lanchester-type  equations 
are  a  gocd  approximation  to  the  mean  course  of  combat.  It  is 
clear  that  ths  casualty  rate  is  equal  to  the  reciprocal  of 
the  expected  time  for  a  force  to  inflict  a  casualty,  when 
the  times  between  casualties  are  exponentially  distributed. 
However,  in  case  the  times  between  casualties  are  no  longer 
exponentially  distributed,  Taylor  used  the  suggestion  made 
by  Bonder  and  3arfoot.  Bonder  [Ref.  3]  and  3arfoot  [Ref.  5] 
suggested  defining  the  Lancnester  attrition-rate  coefficient 
as  the  expected  rate  which  a  single  firer  Kills  enemy 
targets. 

In  the  spirit  of  Bonder  and  Barfoot,  Taylor  provided 
more  rigorous  justification  [Ref.  7]  for  the  Lanchester 
attrition-rate      coefficient      that      does      not    assume    an 
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exponential  distribution  for  times  between  casualties.  He 
considered  the  case  where  the  initial  force  size  of  X  and  Y 
is  large  enough  to  insure  a  negligible  probability  that  the 
battle  is  terminated  before  an  attrition  pattern  could  be 
established-  He    made      no    specific      assumptions   about      the 

distribution  of  times  between  kills,  but  assumed  that  each 
individual  Y  force  kills  X  targets  according  to  an  attrition 
process  in  which  the  times  between  Kills  are  independently 
and  identically  distributed  random  variables.  Thus  in  the 
parlance  of  the  theory  of  stochastic  processes,  he  said  that 
such  an  attrition  process  is  called  a  renewal  process. 
Prior  to  this,  attrition  prediction  nas  been  difficult  due 
to  the  inability  to  predict  casualty  patterns.  By  investi- 
gating the  distribution  of  the  times  between  casualties  via 
renewal  theory,  combat  attrition  analysis  becomes  easier. 
Indeed,  we  know  a  great  deal  about  tne  casualty  pattern  once 
we  know  the  distribution  of  the  interarrival  time  of  kills. 
As  we  consider  more  variables  from  the  renewal  process,  we 
find  more  information  about  tne  casualty  pattern. 
Furthermore,  we  can  track  the  development  of  tne  combat 
attrition  more  precisely  by  studying  the  probabilistic 
distribution  function  of  those  random  variables  that  can  be 
considered    in   the   renewal  process. 

We  assume  here  that  the  reader  is  aware  of  the  concept 
of  a  counting  process.  When  we  think  about  combat  between  X 
and  Y  forces,  we  can  assume  that  time  between  kills  by  a 
single  firer  has  some  distribution  F.  Now  we  are  interested 
in  observing  the  occurence  of  casualty  and  the  number  of 
casualties,  N(t)  ,  that  have  nappened  in  time  interval  (0,t]. 
Furthermore,  when  times  between  kills  are  i.i.d.,  we  call 
this    counting   process  as  renewal    process. 

In  this  chapter,  the  study  of  attrition  analysis  will 
include    the   following    : 
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•  The    r.v.,      N  (t)  ,      total    number      of    casualties    by    tine    t, 
it's   distribution,    and    properties. 

•  The    expected   number   of    Casualties   by    time    t. 

•  Theorems  discussing    casualty    occurence. 

•  The    limiting   behavior   of    the    casualty    occurence. 

•  The   age   and    excessive   life    of    the    casualty   occurence. 

•  Examples  with   specific    casualty    distributions. 

•  Application    of    renewal    process   for   casualty    estimation. 

•  Total      force     attrition      by      pooling      the      single      firer 
renewal    processes. 

B.       REVIEW    OF    RENEWAL   PROCESS 

If  the  sequence  of  nonnegative  random  variables  {X  ,X  , 
...,}  are  i.i.*d.r  then  the  counting  process  [N(t)f  t>0  }  is 
said  to  be  a  renewal  process  (or  often  called  an  ordinary 
renewal  process) .  Thus,  a  renewal  process  is  a  counting 
process  such  that  the  time  until  the  first  event  has  seme 
distribution  F,  the  time  between  the  first  and  the  second 
event  has,  independently  of  the  time  of  the  first  event, 
the  same  distribution  F,  and  so  on.  So  waen  a  event  occurs, 
we    say   that  a   renewal  has   taken   place. 

An  example  of  a  renewal  process,  let  us  suppose  that  we 
have  an  infinite  supply  of  lightbulbs  whose  lifetimes  are 
i.i.d..  Suppose  also  that  we  use  a  single  lightbulb  at  a 
time  and  when  it  fails  we  immediately  replace  it  with  a  new 
one.  Under  these  conditions,  {N(t)  ,  t>0  }  is  a  renewal 
process  wnen  M (t)  represents  the  numner  of  lightbulbs  that 
have  failed  by  time  t.  Accordingly,  if  we  assume  that  this 
failure  is  a  kill  by  a  single  firer  -  X  or  Y  -,  then  killing 
process   also  can   be    considered      as   a   renewal    process.      Thus, 
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renewal     (or   event   or   arrival)         will   be    synonymous    with    kill 

or   casualty   occurence      in   the   rest   of    this      thesis.         He   now 

assume    that    ca s u al t y    occurence   is   a    kind   of   renewal    process. 

If,    in   particular,    the   distribution   of    interarrival   time   of 

kill    is      exponential      with    probability   distribution   function 
x 
e  ,    then  this    renewal    process   turns   out    to    be    a    Poisson 

process    with  rate,       .      [  Ref  -    9] 

Thus  we  may  assume  that  X  casualty  occurence  due  to  a 
single  Y  firer  who  uses  the  same  acquisition  mode  to  engage 
each  X  target  can  be  explained  by  ordinary  renewal  process 
because  we  assume  that  times  between  casualties  are  all 
i.i.d.  However,  we  often  consider  a  counting  process  for 
which  the  first  interarrival  time  has  a  different  distribu- 
tion than  the  remaining  ones.  So,  we  may  think  of  a  Y  firer 
who  at  first  uses  serial  acquisition  mode  to  engage  an  X 
target,  tut  once  he  has  finished  engaging  an  X  target,  then 
he  may  use  parallel  mode  to  engage  the  rest  of  the  targets. 
So,  in  this  case,  we  can  assume  that  interarrival  time  of  X 
casualty  has  different  distribution  between  the  first  and 
the    rest. 

Formally,    let    {X    ,i>1}     be    a    sequence    of    independent   r.v. 

i 
with    X      having    distribution    G,    ana    [X    ,i>2}    having    distribu- 

1  i 

tion    F   which  is    different   from   G.         Then   counting    process 

[N    (t),t>}     is   said    to  be   a    delayed   renewal    process.         So,    we 
D  

can    think      of   using      delayed  renewal      process    for      different 

acquisition   mode,    but   here    in    this    thesis,    we    will    only   deal 
with    the  case    of   ordinary   renewal   process    for   simplicity. 

We  now  define  for  the  renewal  process  a  number  of  asso- 
ciated random  variables  whose  study  is  the  objective  of  this 
thesis-  Table  IV  and  Figure  3. 1  give  the  intuitive  interpre- 
tation  of   the   random    variables    obtained    in    renewal    process. 

Here    we    can   interpret  each    random   variable    as   following: 

•    X      is    the   time    between    the    (n-1)st    and    n-th    X    force   kill 
n 
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TABLE  IV 
Random  Variables  in  Renewal  Process 


r.  v  interpretation 


X        time  between  the  (n- 1)  st  and  n-th  renewal 
n 

i.e.  n-th  interarrival  time 

S        time  up  to  n-th  renewal 
n 

N(t)      total  number  of  renewals  in  (0r  t] 

M(t)      expected  number  of  renewals  in  (0,  t] 
known  as  renewal  function 


Z  (t)      length  of  time  measured  from  the  last 

renewal  until  a  given  time  t  ,  called  as 
age  or  backward  recurrence  time 

Y  (t)      length  of  time  measured  from  time 

t  until  the  next  renewal,  called  as 
excess  life  or  forward  recurrence  time 

S(t)     length  of  renewal  time  at  time  t, 
called  as  spread 


by  a  single  Y  firer. 

•  S   is  the  time  up  to  n-th  kill. 

n 

•  N  (t)  is  the  total  number  of  kills  in  (0,t]. 
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< s  (t)     > 


< X >< X X- 

1  2 


<— Z(t)— > 


N(t) 


<— Y(t)~  > 


time 
•-x 

'N(t)  +1 


Figure  3.1         Renewal   Process. 

•  M  (t)    is   average    number    of    kills    in    (0,tj. 

•  Z  (t)    is   the    time   between   the    last    kill    and  a    given    time    t 

•  Y  (t)     is   the    length   of   time   from   time    t    up   to    the 

next   kill. 

•  S  (t)  is  tae  length  of  lifetime  of  X  force  soldier  at 

time  t. 

In    this    thesis,       we    will   only   deal    with   X   force    casualty 
by    Y    firer    to   avoid    the   confusion   in    notation. 

C.       DISTRIBUTION    OF    TOTAL    NUMBER    OF    CASUALTIES 

Suppose  that  for  an  ordinary  renewal  process  a  single  Y 
firer  uses  the  parallel  acquisition  mode  to  engage  passive  X 
targets  from  tae  beginning  of  combat,  when  we  assume  that 
targets  are  alrsady  acquired.  Then,  the  distribution  of 
N  (t)  ,  the  total  number  of  X  force  casualties  by  a  single  Y 
firer,  can  ba  obtained  at  least  in  tneory  by  first  noting 
the  relationship  that  trie  number  of  kills  bv  time  t  is 
<2£§^ter  than  or  e<qual  to  n  if  and  onlv  if  the  n-th  kill 
occurs  before  or  at  time  t. 
Formally 
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N (t) >n    <==>    S    <t  (3. 1) 

n 

From    3.1,    we  attain 

?{N(t)=n}    =    P{N(t)>n}    -   P{N(t)>n+1} 

=   P  (S    <t}    -    2{S         <t} 
n  n+1 

Now,    since    the    random   variable   X    •s#(i>1),    are   independent 

i  n 

and    have   a    common   distribution    F,    it    follows   that    S    =  IEL  X 

n      i=1    i 
is   distributed    as   F    ,    which    is   n-fold    convolution    of    F    with 

n 
itself   [Ref.    10]-         Here,       F   is    tne    probability    distribution 

of    time    retween    casualties.       Therefore   we    obtain 

P[N(t)=n}    =    F     (t)     -   F         (t)     with    F     (t)     =    1 

n  n+1  0 

and 

P{N(t)  <n}    =   ?{S   >t)    =    1    -    F     (t) 

n  n  (3.2) 

Therefore    the      probability    distribution      of    N (t)  ,         which    is 

total      number      of   casualties      by      single      Y    firer,         can      be 

obtained   explicitly    fcr   all    n. 

The    simplest   special   case    of   equation    3.2   is    obtained   by 

taking    the    renewal    process    to    be    a    Poisson    process. 

Then    S      has   the      special      Erlang      distribution   with    n    stages 

n 
because    each  F       is    exponentially   distributed.      Also,         it   is 

known  from      the    previous   argument      that    N(t)       has      a    Poisson 

distribution  with   mean        t. 


L.       EXPECTED   30*flBER    OF    CASUALTIES 

1 .      Renewal    Func  tion 

Expected  number  of  casualties,  E[N(t) ],  which  is 
defined  to  be  the  mean  number  of  casualties  in  the  interval 
(0,t]#    is   known    as    the   mean    value   or   renewal   function,    M(t). 

From    equation   3.1,    the   expected   number    of    casualties   in   time 
t  can    be    obtained  as    following: 
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3<t) 


=    E[N(t)  J   «  J     k   ?[N(t)=k} 

=ZL   21     P{N(t)=k} 
n=  1  k=n 

fJO 

=  IP(N  {t)>n} 
n=1 

=  2_   P[S    <t} 
n=1        n 

bo 

=  X  F    (t) 
n=1    n 


(3.3) 


Eut   here  it    is   generally    very    difficult    to   find    M (t) 
directly  from   the   distribution   F.       Thus,       we    will    employ   the 
Laplace    Stielt  zes   Transform     (L.S.T.  )         [Ref.    9]    to    find  .more 
conveniently      an   expression      for    H(t) .  We    generally      need 

three    steps   to    find    M  (t)  . 

Ste£    1: 

From   the  definition    of    L.S.T. ,         we   know    tnat    L.S.T. 

of    a  (t)    is 

*  f        -st 

M  (s)    =        e  dM(t) 

0 


and   also 


F(s)     = 


-st 
e  dF(t) 


Step    2: 


Then  from   equation    3.3 


oo 


H(S)     =-      F     (s) 

n=  1      n 

'^  *         n 

=  ^  :f(s)  ] 

n=1 


for    s>0 


The  above  is  a  sum  of  geometric  series  and  can  be  written  as 
follows: 


*       F(s) 

a  (s)  =  

1-F(s) 


(3.4) 
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or   equivalently 

*                 M(S) 
F(s)     = -— 

1+M(S) 

Step    3: 

Since  the  L.S.T.  of  a  function  uniquely  determines 
the  function,  the  renewal  function  determines  the  distribu- 
tion of  the  kill  time  and  therefore  determines  the  prob- 
ability law  of  the  casualty  occurence.  Specifically  there  is 
a  one-to-one  correspondence  between  interarrival  time 
distribution  F,  and  renewal  function,  II  ( t)  .  So  if  we  have 
the  distribution  F,  then  we  can  find  d (s)  from  the  equation 
3-4.  But  to  find  tf(t),we  have  to  transform  H  (s)  to  M  (t)  by 
inverse   application    of   L.S.  T, 

Let  as  try  the  example.  Remember  that  this  example 
will    be    applied    all    the    way    through   this    chapter. 

Example   3. 1 

Suppose    that    each   interarrival   time    to    kill    X    ,(i>1)    has 

i 
exponential   distribution  witn      .    Then    L.S.T.    of    F    is 

(*    -st 
F(s)     =  e         dF(t) 

f*     -st         -H 

=    I       e  /  e        dt 

-  0 

-{/  +  s)t 

=        Ae  dt 

0    ' 

a        />0  -  (A+s)  t  — 

=      (A+s)  e  dt 

A+s        0 


A 
A+s 


if   s>-'A 

Then    from   equation    3.4 
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*  F(s)  A/(A+s)  >. 

M  (s)      = ■     = = 

1-F(s) 

where  1/s  was  an  L.S.I.  of  F (t) ,  which  had  value  of  t.3  So 
M(t)  =  ht  +  c  (some  constant)  ,  bat  we  know  that  8  (0)  =  0, 
because  there  would  be  no  casualty  at  time  0.  So,  finally 
M(t)     =  At. 

Let  us  be  more  specific  about  this  example.  Suppose 
that  interarrivai  times  of  casualties  are  i.i.d.  and  expo- 
nentially  distributed   with    mean    30    minutes.    Then 

a)  Wnat      is    the    distribution   of      total    X   casualties      at   5 
hours   later    ? 

b)  What    is   the   probability    that    3rd    X   kill    was    occured   at 
time   t    =  4    ? 

c)  What   is  the    prctability   that    at    time    t   =    2,       no    X    casu- 
alty  occured    ? 

To  answer  these  questions  we  know  that  M (t)  =  2t  from  above 
example  because  =  2,  then  number  of  X  casualties  in  time  t 
has    Poisscn   distribution    with    mean    2t. 


3Let    F(t)    =    t   for   t>0.     Then    L.S.T.    of    F(t)    is 

*                f°°    -st                   C*     -st                        1         -st  ^ 
F  (t)    =      \      e        d?  (t)    =        e        1dt   = e 

0  0  s  -0 

=    1/s   for   s>0 
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Then 


-10  n 

e  (10) 

a)     P{N(5)  =n}    = -  — 

n    ! 


3      -8 
8      e 
b)    P{N(4)=3}    ■ =    .023  6        and 

-J        * 


0      -4 
4      e 

c)     ?{N(2)=0}    = =    .0183 

0    J 


2.      Benewal    Equation 

Now  an  integral  equation  for  M(t)  may  be  obtained  by 
conditioning  on  the  time  of  the  first  casualty  occurence. 
Upon    doing    so,     we  obtain 


( 


JO 


H(t)=  E[N(t)  |X^   x]   dF(x)  (3.5) 

ho  w  e  v  er 

<  0  if   x>t 

E[  N(t)  |X    =    x  ]    =      i 

1  v.  1    *   M  (t-x)  if    x<t  (3.6) 

for  if  the  first  casualty  occurs  at  time  x,  x<t,  then  from 
this  point  en  the  process  starts  over  again,  and  thus  the 
expected  number  of  casualties  in  (o,t]  is  just  1  plus  the 
expected  number  of  casualties  in  time  t-x  from  the  beginning 
of  an  equivalent  renewal  process.  Putting  equation  3.6  in 
3.5    yields 

•t 
M(t)     =    \       :  1  +  M(t-x)  ]dF(x) 

t 
=    F  (t)     +    ■         M  (t-x)  dF  (x)  (3.  7) 
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where  equation  3.7  is  known  as  the  renewal  equation  and  F  is 
a  known  function  and  H  is  an  unknown  function  to  ce  deter- 
mined as  a  solution  to  the  integral  equation  3.7.  Ross 
provided  the  solution  for  equation  3.7  [Ref-  11]  as 
following  and  the  solution  is  accepted  true  for  every  case. 

Expected  number  of  casualties  up  to  time  t  is 


M(t)  =  F(t)  ♦     F(t-x)dH(x)  (3.8) 


where  H  (x)  =  ZL   F  (x) 
n=1   n 

So   it    allows  to    find   M(t)    using    only    the    known   distribution 
function   F,    which   is    the   distribution    of    time   between   casualt 

ies    ,    thus    avoiding    finding    individual   F    . 

n 

3.      Renewal    Density 

Now  let's  define  the  renewal  density  as  m(t)  = 
dil{t)/dt.  Remember  that  M(t)  is  the  expected  number  of  casu- 
alties up  to  time  t.  So  a(t)  may  be  regarded  as  the  expected 
number  of  casualties  per  unit  time;  more  precisely,  the 
average  number  of  casualties  during  the  time  interval  from  t 
to    t+h   is 

't+h 

M(t+h)     -    M(t)     =  \  m(x)dx  (3.9) 

Differentiating  the  renewal  equation,  we  obtain  the  renewal 
equation   for   the    density   m(t) 

.,.,.«,.(>.,.,„*  «„„. 

So, this  knowledge  of  the  density  f  of  interarrival  time  of 
casualties  is  sufficient  to  determine  the  density  m  (x)  of 
casualty   occurence. 
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E.       LIMIT    BEH&VIOOB    CF    CASUALTY    OCCURENCE 

3e  know  that  total  number  of  casualties,  N  (t)  ,  go  to 
infinite  as  time  goes  to  infinite,  when  we  assume  that  force 
size    is   infinite-  Eut  it    would    be    useful   to      know    at   what 

rate  N(t)  approaches  infinity.  That  is,  we  would  like  to  be 
able  to  say  something  about  limit  of  N (t) /t  as  time  goes  to 
infinite. 

From   the  definition   of    r.v.    we   know    that 

S  <    t    <    S 

N  (t)  M  (t)  +  1 

wnere   S  must    be   the  time    of    the    casualty   after    time   t. 

N  (t)  +1 
By   refer   to   [Ref.    10]  we   know   that 

N  (t)  1 

— >      as   t   -->  Oo  (3.11) 

t  u 

where   u    =  E[ X        =     l        [  1    -    F(x)  ]dx    ,    wnich    is    the 
i  0 

expected    time    between   casualties. 

Hera  the  function  1/u  is  often  called  as  the  rate  of  the 
renewal  process.  Thus  the  average  number  of  casualties  per 
unit  time  converges  to  1/u  .  Then,  how  anout  the  expected 
average  number  of  casualties  per  unit  time?  Is  it  true  that 
H(t)/t  also  converges  to  1/u  as  t  -->  0&  ?  This  result,  known 
as  the  elementary  renewal  theorem,  will  be  stated  without 
proof    [Ref-    11  ]. 

Elementary  Renewal  Theorem 

M  (t)  1 

> as   t-->    \X>  (3.12) 

t  u 


Thus  whatever  the  distribution  F   is,   the  average  number  of 
casualties  up  to  time  t  is  approximately  (for  large  t) , 
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M(t)  = (3.  13) 

u 

Also,  the  average  number  of  casualties  in  the  inter  yal 
(%.*    t+h)  for  h>0  is  approximated  for  large  t  as  following, 

h  A 

M(t+h)  -  M(t)  = as  t  — >  PQ  (3.14) 

u 

for   any    distribution    F   with    mean    a. 

Here  equation  (3.14)  is  known  as  Black,  well f  s  theorem  when  F 
is   not   lattice. 

From  the  example  3.1,  where  casualty  occurence  time  is 
exponentially  distributed  with  mean  30  minutes,  the  average 
number  of  X  casualties  up  to  time  t  is  approximated  by  t/u, 
which  is  2t.  Then  after  2  hours  of  combat,  we  can  assume 
that  the  average  number  of  X  casualties  by  single  '£  firer 
will    be    2t,    which   is    4   casualties. 

Another  limiting  result  which  may  prove  to  be  useful  in 
renewal  process  also  concerns  the  r.v.  N  (t)  .  It  is  shown 
[Ref.  12]  that  N  (t)  has  an  asymptotically  normal  distribu- 
tion   with   mean    t/u    and  variance   t  X^3  .       Ihus 

_  ,■     y   -    t/u   N 

lim      F{N(t)<y}    =g>    .--» -  (3.15) 

t->E*  ~    i      7.  3 

\\G  t/u 

Suppose  of  the  same  example  3. 1  that  casualty  occurence  mean 
time  is  30  minutes-  Then  u  =  0. 5  =  E[  T  ]  =  1 //K  .  So  A  =  2. 
Also  Var£T]  =  1/A-  =  1/4  =  0.25  hr  ,  which  is  900  mint  here. 
Let's  assume  that  we  had  a  10  hours  of  combat  action.  Then 
expected  number  of  X  casualties  is  E[  N  ( 1 0)  ]  =  t/u  =  10/0.5  = 
20    and   var[N(10)  ]  =    t  xVu*  =    10x0-25/(0.5)=    20. 

Suppose  of  the  question  that  how  many  X  soldiers  will  be 
required  in  order  to  with  probability  of  95%  that  at  least 
more  than  one  X  soldiers  will  survive  after  10  hours  of 
combat   ?    This    question   may    be    answered    as    following, 
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_/  X    -    20    \  ->-  '    x    -    20    \ 

0.95    <    P{N(10)<x}     s  cb  —  ~ -— —  '       =>•- --- 

-1     i  10  /  —         4.4  72 


Now  from  the  normal  table,  0.95  percentile  of  normal  distri- 
bution is  1.645.  So  1.645<  (x  -  20) /4. 472,  where  the  answer 
turns  out  to  be  x>27.36  .  So  we  can  say  that  X  force  needs 
at  least  28  soldiers  as  an  initial  force  size  to  te  95*  sure 
that  at  least  more  than  one  X  soldier  can  survive  after  10 
hours    cf   combat. 

There      also    exist     a     alternative   normal     approximation, 
when    we    assume    that    for   integer   n, 

n 
P{N(t)<n}    =    P  [S    >t}      where    S      =       7Z     X 

n  n  i=1      i 

from    equation  3.  1. 

_  /   t  -  nil   \ 

Then    P{S    >t}    =     1    -    P  {S    <t}     =    1    -    £.  j -rjr- '  (3.16) 

n  n  ry-'i  n 

So   now   we   assume    that  S      has   normal    distribution    with    mean 

n 
nu   and   variance    n  y  .      If      we   apply   equation    3.16    to    previous 

example, 


-_  I  10    -    0.5x    \ 
.  95<P{N  (10)  <x}     =    1    -  $£_ ZZT-! 


0.  5  j  x 

—   /  10   -    0.5x    ,  10    -    0.5x 

Then  <£>  ; -____;<    q.05      or      =rr-      <   -     1.645 

"" '  •  0.  5  'HT   /  0.  5  J  x 


If  we  solve  for  x,  it  turns  out  to  be  x>28.83.   So  answer  is 
almost  same  when  we  assume  some  rounding  errors. 

F.   AGE  AND  EXCESS  LIFE  DISTRIBUTION  OF  CASUALTY  OCCURENCE 

As  we  mentioned  earlier  in  this  chapter,  we  can  also 
consider  such  random  variables  as  Y  ( t)  ,  which  is  the  time 
from  certain  time  t  until  next  casualty  occurence, 
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which  means  residual  lifetime  of  a  soldier  in  comhat  from 
certain  time  t  if  he  is  still  alive,  and  Z  (t)  ,  which  is  the 
time  from  t  since  the  last  casualty  occurence  (see  Figure 
3.2). 


Figure   3.2        Age    and   Excess    Life   of   Casualty  Occurence. 


This    means 


y  (t)    =  s  -   t 

n  (t )  +  1 


Z  (t)    =   t  -  s 


N  (t) 


where  Y  (t)  is  called  as  the  excess  life  or  forward  recur- 
rence tiae  and  Z  (t)  is  called  as  a£e  or  oackwar a  recurrence 
time.   Also 

S(t)  =  Z(t)  +  Y(t) 

=  S        -t+t-5 

N(t)  +  1  N  (t) 

=  S        -  S 

N(t)  «•  1     N(t) 


N  (t)  +1 


& 


|tt, 


X   -    27   X 
i=1    i      i=1   i 


=  X 


N(t)+1 


(3.17) 
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Here    S(t)    is  often    called   as   spread   and    S(t)     is   also    r.v. 

because    X  is    r.v,    too - 

N(t)  +  1 

Suppose    of    the  Poisscn    process    that    X.    is    exponential, 

then    E[S(t)  ]  =    E[  I  (t)  ]   +    E[  Z  (t)  ]    =    1/A    +    V,\  =    2/\    =    2E[X.]. 

So,    S  (t)    can  be    assumed   as    the   length    of    the    lifetime    of    the 

X  force    soldier    which  is   in    combat    at    time    t. 


IV.    THE    SU PER POSIT I  ON    OF    CASUALTY    OCCURENCES 

A.       INTRODUCTION 

In  the  previous  chapter,  we  have  considered  in  detail 
the  theory  of  casualty  occurence  through  the  ordinary 
renewal  process  of  a  single  Y  force  soldier.  We  deal  in  this 
chapter  more  briefly  with  the  superposition  of  several 
renewal  processes,  which  is  total  X  force  casualty,  occurence 
behaviour.  When  we  assume  multiple  Y  force  soldiers  who 
shoot  at  X  force  soldiers  independently  of  each  other,  we 
may  assume  that  all  these  Y  soldiers  kill  X  soldiers  by 
independent  casualty  processes.  So,  we  can  say  that  there 
occur  multiple  independent  ordinary  renewal  processes, 
simultaneous  in    time. 


process    1    | 


process    2   \ 


process    3   J 


pooled 
output 


x x x 

x ' x ' x ' 

x— ' ' * x ' * x ' 

II  II  III  II 

II  II  III  II 

It  It  III  It 

x--x x x x x x x X 


Figure   4.1        Superposition   of   Casualty  Occurences. 
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Suppose  that  we  have  '  p'  independent  casualty  processes 
in  operation  simultaneously ,  all  with  tne  same  P.p.F.  of 
casualty  occurence  time.  Consider  tae  sequence  of  casualty 
occurences  formed  by  pooling  the  individual  processes.  Fig 
4. 1    illustrates    the    special    case,    when    p   =    3.    Formally,    we 

can    say   that  if    {N    (t),    i>0},      for   i   =    1,2,      ...       ,      k,      are 

inde- 

i  k        pendent        renewal 

processes   and  N     (t)  =        N    (t)  ,    then    £N     (t),t>0} 

p  i=1    i  p     is 

called   pooled  process   which    is    total    X    force    casualty    occu- 

rence      behaviour      by        k      Y      force 

soldiers. 


B.       SCME   GENERAL   PROPERTIES 

1  -      Normal    A  gprq  ximation    of    Pooled   Output 

From  casualty      occurence    processes      in   the      previous 
chapter,    we    assumed    that   there   exist    a    mean  casualty    occure- 

2 

nee    time.    So  let   u      and   (J7       be   the      mean      casualty    occurence 

i  '  i 

time    and      ccrrespondina      variance   by    the    i-th    Y   force    soldier 

•     2 
and    let    u(p)    and     rj-  (pj    be      the   mean      casualty      occurence    time 

and   corresponding   variance    of    the   total    casualty    occurence   by 

k   of    Y   force  soldiers.    Then    the    question   is   how    are    these    two 

values   related    each    other   ? 

From   previous  study,    we    know    that    1)     N     (t)    is    approxima- 

i 
2  3 

tely    normal    with   mean  t/u      and   variance,    (J"    t/u      .       Therefore 

i  i  i 

we    can    say    that    N    (t)    is   also   approximately   normal    with   mean 

k  ?  k  2  3 

TL     t/ui   and  variance     !_      CT~   t/u      .       2)     Also    assume    now    that 
i=1  i=1      i  i 

either      the   pooled    process    is    a      renewal    process    or    at    least 

the    analogous  result    holds,     tnen    approximately   N     ( t)     is   norm- 

2  3  ? 

al    with    mean  t/u  (p)     and    variance      (J  (p)  t/    u    ( p)  . 

So,    from    1)    and    2)  ,     we   must    have 

t  k         1 

=   t    r      

u  (p)  i=1      u  (4.  1) 
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k      1         -1 

then    u(p)    =    ( yi )  which   is    the    harmonic   mean    of    the  u    , 

i=  1    a  i 

i 

2  2 

(P)                     k           i 
and t   =    t     F_     t 

3  i=1  3 
u     (p)                               u. 

l 

2  ^  2 

2  3  k      (Ti  k      1        -3      k     Qi 

then     j-(f)     =   u     (p)       r =    (Sl~-    )  £: --     (4.2) 

1  =  1  3  i  =  1    u  1=1         3 

u  i  u 

i  i 

Also    if    we    assume   that   u      =    u    and  <]~~    =        for    all    i    =    1,2,..,k 

i  i 

which    means   all    casualty   occurences   are    independent   and   identi- 
cal,   then    mean    casualty    occurence    time   and    corresponding    vari- 
ance   of   the   total  casualty      ioehaviour      by    k    of   '/    force    are  as 
follows; 


1 

-1 

u 

u 

(F) 

— 

( 

u 

k 

) 

— 

k 

2  2 

2  1-3    0"  (T  r     2 

rp    (p)    =    (  -    k    )  —  k    =—-=(-    )  (4.3) 

u  3  2  k. 

u  k 


Also  total    number  of   expected  casualties,    E[  N    (t)  ]r    can 
p 

be    obtained   here   by    pooling    individual   casualty    occurences 
as    following; 

k 
E[  N    (t)  ]   =     H      E£  N    (t)  ]    =    k    E[  N    (t)  ]    =    k    M  (t)  (4.4) 

p  i=1  i  i 

where    H  (t)    is    expected   number    of    casualty    from 
single    fixer,     which    is    t/u   for   large    t. 
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let's  look  at  an  example  that  deals  with  pooling  of 
individual   renewal    processes. 

Example    4. 1 

Suppose  of  Y  force  artillery  shells  from  several 
sources  are  bombording  the  same  X  target.  Each  source  hurls 
shells  at  the  target  at  a  rate  of  40  shells  per  hour.  Assume 
that  the  interarrival  times  of  shells  from  each  single 
source  are  uniformly  distributed  over  an  interval  (0,  a]. 
Then  the  question  is  1)  What  is  the  probability  that  less 
than  825  shells  will  be  hurled  at  the  target  from  a  single 
source  during  a  20  hour  period  and  2)  If  there  are  5  sources 
from  which  such  attacks  are  launcned,  what  is  the  prob- 
ability that  more  than  4125  (  which  is  5  times  825  )  shells 
are    hurled    at    the   target   in    a    20    hour    period    ? 

To  answer  the  first  question,  let  A  =  40  shells/hour 
=  2/3  shells/mi n.  Then  u  =  E[  X  ]=  1/-\  =  3/2  min,  when  we 
assume  pcisson  process.  But  interarrival  times  are  uniformly 
distributed  over  an  interval  (0,a],  So,  u  =  E[X]  =  (0  + 
a)/2=  a/2,  where  a  =  3  minutes.  Then  Var[X]  =  (a  -  Q)~/12  = 
3/4  min~  .  Now  let  N  (t)  =  number  of  shells  hurled  in  (0,  t] 
from  a  single  source.  Since  N  ( t)  is  normal  with  mean  t/u  and 
variance  tr/u3  from  equation  3.15,  for  t  =  20  hours  =  1200 
minutes,    the  probability    turns    out   to    be   as    following; 

1200 
825    —    —  —  ——— 

^r  i                       3/2      ]           -^  /  825    -    800    \ 
P  (N(1200)<325}    =  Ct   '-——-— 1       a  a    I 

~~\  I  120  0    -    3/4  —[         16.33 

\~~~(!/2f~~ 
=  £(1.53)     =    0.9370 
and    for   second    question, 

Let    N.(t)    =    number    of   shells   hurled    in    (0    ,    t  ]  by    i-th    source 
where  i    =   1,2,     ...    ,    5    .    and 

5 

B  (t)  «  I  J  (t)  =  total  number  of  shells  hurled  in  (0,  t] 
p     i=1   i 
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by    5    sources. 

Then    {N     (t)  ,   t>0}    is    the    pooled   process    for   which    mean 

p 
inter arrival  time   is, 

u  3/2  3 

u  (p)    =   -   =      =    —    minutes. 

k  5  10 

and   corresponding   variance      is 

2  -  2  3  2 

(TCP)    =     (  Vk   >       =  ----    mm 

Since    N     (t)    is    normal  with    mean   t/u  (p)    and    variance  y  (p)  t/(  i3  (p) 
for    t=   1200    minutes, 

1200 

4  125 , 

-T-  I  3/10   \  -^  /    12  5        \ 

?{N     ( 1200)  <4  125}     =     S     -■ r-  =9 ) 

p  —  ,     1200    3/100  -^A    36.51 


t 


(3/10) 


■(3.42)     =    .9997 


So,     P  {N     (1200)>4125}     =    0.0003 


2-      Poisson    Process   of    Pooled   Out  put 

Now    let's    assume    that     individual     casualty    processes     are 

probabilistically    identical    witn   common    underlying    C.D.F.    of 

casualty      occurence      time      F  (x) .    Let    F     (x)       be      the    C.D.F   of 

p 
casualty    occurence    time    of    pooled   process.    Tnen    the    question 

is    how   is    F    (x)       related      to   F  (x)     ?        This      question      can   be 

p 
answered   by   using   the   idea    of    excess    lire.    Let   Y     (t)       be    the 

i 
excess   life    of    an   X    force   soldier    and    Y     (t)    be        excess   life 

of    total   X    force.    Then   we   know    that    Y    (t)     =   min[  Y     ft),    Y    (t) 

....  ,  Y  (t)  ]  (see  Fig  5.1). 

k 
Then 
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Figure    4.2        Superposition   of    Excess   Life. 


?{Y    (t)>y)    =    ?{Y.  <t)>y} 

P  i 


(4.5) 


From    previous   study,    we    know   that    P.D.F   of    E.D.    for    individual 

renewal    processes  is    (1/u)  [  1    -   F  (x)  ]    and   then    we    can    say    that 

1 

[1  -  F  (x)  ]  is  P.D.F  of    E.D-  for  the  pooled  casualty 

u(P)  P 

process.    Thus, 
P£Y(t)>y}    = 


,-CO     1 

-    [  1    -   F  (x)  ]dx      then, 
y     u 


f          1 
P{Y    (t)>y}    =     \  [  1    -    F     (x)  ]dx 

P  'y      u(p)  p 


)y      u(p) 

1     /  k 

=  [     -     \      [  1    -    F  (x)  jdx]      from   equation    4.5 
u      I  y 


If    we    differentiate    with   respect    to    y 
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1  1  "*>  k-1    1 

[  1    -   F     (y )  ]    =  k[    -    ,     [  1    -    F  (x)  ]dx  ]       (-   -    )  [  1    -   F  (y)  ] 

u(p)  p  u    ' y  u 

Here    we    know   that  u(p)    =   k/u,    so    u's    are   all    cancelled    cut. 
Then, 

C.D.F.    of   pooled   output   is 

^     r  k-1 


p  u        y 

For   example,    think   of    the   exponential    distribution   of 

-  X  x 
casualty    occurence.    Then   F(x)    =    1    -    e  ,    and   u    =    1/\    . 

Then    C.D.F    of    pooled   casualty    occurence    is 

^    -ax        k-1  -Ay 

F     (y)    =    1  -    [       \     [e        ]dx]        e 

p  y 

->.y   k-1   ->y 
=    1   -     (e         )         e 

-Xky 
=    1  -    e  for   y>0 

which    follows    that    pooled   output   is    a    Poisson    process    with 

rate    k      where   k    is    the   number    of    1   force   combatants. 

Here,  probably  the  most  interesting  properties  of 
the  pooled  output  refer  to  the  'local'  behaviour  where  k  is 
large.  Khintchine  [  Ref .  13]  has  proposed  that  in  the  limit 
the  numbers  of  renewals  in  non-overlapping  intervals  follow 
independent  Poisson  distributions,  thus  showing  that  in  the 
limit   the      pooled  output     is  a      Poisson    process.  Also   his 

proof  does  not  require  each  interarrivai  time  to  be  identi- 
cally distributed.  In  attrition  analysis  in  combat,  it 
appears  to  be  the  same  way,  which  means  that  if  a  large 
number  of  independent  (and  identical)  casualty  occurences 
are  pooled,  then  the  total  force  casualty  occurence  is 
approximately   a    Poisson   process      [Ref.    9]   [Ref.     13].         Thus, 
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we  can  say  that  pooled  casualty  occurence  by  total  Y  comba- 
tants can  be  approximated  b  v  ggisson  process  if  we  assume  a 
larje    d "^ gb § r  of    Y   forces. 

Thus,      finally,         from    the    nasic    theorem      of    Poisson 
process,    we  know   that  if    {N (t) ,t>0}    is    Poisson  process    with 

rate      >0,    then    the    interarrival   times    {X    ,i>1}    are   i.i.d- 

i  -, 
exponential     candom      variable      with    rate    A    .         Thus,         from 

previous  results,      total  X    casualty      occurences   is   a   Poisson 

££2cess   with     rate   Jc /V  ,        where   k   is      the   number    of      Y    force 

combatants,      and  interarrival     times   of    X   casualties      can   be 

considered  to  be  exponentially  distributed  with  parameter  k 
A  " 
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V.     COMPARISON    OF    RENEWAL    PROCESS    WITH    COMPUTER    SIMULATION 

In  this  chapter  we  will  Iook  at  the  results  of  seme 
typical  combat  attrition  process  (  Y  firers  with  different 
cumber  of  passive  X  targets)  obtained  through  the  theoret- 
ical renewal  process  model  which  is  developed  in  this  thesis 
and    computer   simulation    using    random    number   generator. 

As  we  mentioned  earlier,  total  number  of  expected  casu- 
alties by  n  1  firers  turned  out  to  be  n  M(t)  for  large  t 
from  the  renewal  theory,  where  we  assumed  that  all  casualty 
occurences  were  independent  and  identical.  But  it  is  clear 
that  these  expected  casualties  can  not  be  bigger  than  the 
number  of  X  targets.  Now  we  will  compare  this  result  with 
simulated  output  to  see  whether  we  can  find  any  interesting 
facts  such  as  how  many  X  targets  will  be  most  appropriate  to 
apply  renewal  theory  to  combat  attrition  ?  or  what  force 
ratio    will    be   most    applicable    to    renewal   process    ?,    etc. 

In  order  to  facilitate  the  simulations,  the  assumptions 
were  made  that  n  firers  shoot  independently  at  m  targets, 
and  that  the  times  to  kill  the  targets  are  uniformly 
distributed   between    [0,1].  When   a    target   is      killed,      all 

firers  shooting  at  it  are  assumed  to  select  a  new  target  at 
random   from   those  stiil   surviving. 

The  procedure  of  the  computer  simulation  is  as  follows 
(see    Appendix   A)  ; 

1.  We  generate  the  initial  kill  time  for  each  firer  by 
using  pseudo  random  number  generator,  where  we  assume 
that  every  firer's  kill  time  is  uniformly  distributed 
between    [ 0, 1  ]. 

2.  Each   target    is  indexed    by    integer    number. 

3.  We  need  to  decide  which  target  does  the  firer  snoot 
at   ?      The    target   selection   rule   is    such  that    first    we 
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generate  a  random  uniform  number  between  [0,1],   then 

we  take  the  integer  value  of  (random  number  times  the 

number  of  targets  )   plus   one,   which  designates  the 

index  of  target   that  is  assigned  to   firer  1.    Same 

procedures   are  repeated   until   all   the  firers   are 

assigned  a  target,   with  some  targets  possibly  hcving 

more  than  one  firer.    Thus  at  first  time  step  random 

number   for   selecting   the  target   is   generated   n 

times,   because   we  have  to   assign  a  target   to  each 

firer.    These  targets  are   killed  eventually  as  time 

proceeds. 

Now,       the      target   that      will  be   killed      first    is      the    target 

which   is   assigned      to   the   firer    who   has    the      minimum   time   to 

kill  (   =    Tmin)    among    the   firers,      where   this   time    to   kill   was 

already    generated.    So   we   now   know    which    target   is    killed. 

After  the  first  target  is  killed,  the  number  of  targets 
killed  is  registered  as  one  and  the  first  time  step  ends. 
Now,  in  the  second  time  step,  tne  new  kill  times  between 
[0,1]  are  generated  again  for  all  firers  who  were  shooting 
at    the      first    target.  The   firers      who   were      still   in      the 

process  of  killing  a  target  at  tne  end  of  the  first  time 
step  continue  to  engage  the  same  target  in  the  second  time 
step,  where  time  to  kill  in  the  second  step  is  old  kill  time 
(   generated   value   at    first    time    step    )     less   Tmin. 

Cnce  the  second  target  is  killed,  the  second  time  step 
ends  and  the  number  of  targets  killed  is  now  registered  as 
two.  Now  we  proceed  to  the  third  time  step,  where  same 
target  selection  rule  and  firing  rule  are  applied  until  all 
the    targets   are    killed. 

After  all  the  targets  are  killed,  each  time  step  length 
(  which  is  the  same  as  the  time  to  kill  the  target  in  each 
time  step  )  is  stored,  and  this  is  one  replication  (see 
Appendix  E) .  To  obtain  the  expected  value,  40  replications 
of    the    simulation      were    conducted    and   6    data    sets 


46 


(11=3,  111=6),  (0  =  3,111=10),  (l=3,  m=14)  ,  ana  (n=1G,  m  =  6)  ,  (n=10, 
m=10),  (n=10,  m=  14)  ,  were  examined.  After  40  replications  , 
the  data  were  evaluated  to  find  the  expected  number  of  kills 
through  order  statistics.  This  expected  numoer  of  kills  was 
plotted  against  time.  The  results  for  n=3  and  n= 1 0  are 
shown    in   Fig  5.1. 

The  theoretical  results  can  be  obtained  through  two 
methods,  whi^h  are  discrete  approximation  and  simulation. 
In  this  thesis,  we  used  discrete  approximation  based  on  the 
renewal  function  modelled  at  equation  3.7.  It's  formula  is 
as    follows; 


M(t)     =     21  [     1    +    M  (t-j)     ]    P. 
3  =  1  3 

where    P      represents    the    probability    distribution 
of  kill   time. 


Thus,      expected      number   of    kills      obtained    from      the    renewal 
process    model      using    discrete    approximation    (see      appendix   C 
for    program   and    data,      where   left      column    refers    to    the   kill 
time    and  right    column   refers    to   the    expected    number   of    kills) 
is         alsc    plotted   against    time    in   Fig    5.1. 

It  was  found  through  both  cases  that  the  expected  number 
of  kills  by  the  renewal  process  model  may  only  be  justified 
through  the  assumption  of  an  infinite  numner  of  targets 
because  as  m  gets  bigger,  the  expected  kill  value  approaches 
the  theoretical  result.  Also  it  was  discovered  tnat  at  the 
beginning  of  combat,  the  attrition  process  almost  followed 
the    theoretical    renewal    process    model. 
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THEORETICAL  VS  SIMULATED  RESULTS 


3  RRERS 


10  FWERS 


TIME 


Figure   5.1        Expected  Kills    From   Theoretical   And    Simulated. 
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VI.    CONCLUSIONS 

A,       SUMMARY 

Since  Lanchester  proposed  his  model,  many  studies  have 
tried  to  model  the  combat  attrition  process  mathematically. 
Recently,  the  main  focus  has  been  on  calculation  of  numer- 
ical values  far  Lanchester  attrition-rate  coefficient.  This 
numerical  determination  stands  at  tne  heart  of  casualty 
assessment  in  such  models.  There  are  three  reasons  for  the 
importance  of  this  attrition-rate  coefficient  methodology. 
First,  Lanche ster-t ype  models  are  used  in  various  U.S. Army 
and  D.O.L.  planning  activities  more  widely  than  ever  before. 
Second,  Lanchester  attrition-rate  coefficient  is  a  basic 
element  of  any  Lanchester-type  combat  model,  and  that 
attritioD-rate  coefficient  reflects  tne  effective  applica- 
tion of  firepower.  finally,  significant  new  developments 
have  occured  in  methodology  for  developing  more  tactically 
realistic  Lanchester  attrition-rate  coefficient  and  these 
important  results  have  not  been  accessible  to  a  very  wide 
audience.  In  particular,  a  new  approach  for  developing  more 
realistic  and  mathematical  Lanchester  attrition-rate  coeffi- 
cient, that  of  computing  combat  attrition  by  using  a  renewal 
process,  was  proposed  by  Taylor  recently.  This  thesis  began 
with  the  assumption  that  the  combat  attrition  process  was  a 
kind  of  renewal  process,  and  that  times  between  casualties 
were  i.i.d-  Thus,  considering  the  problem  from  the  stand- 
point cf  renewal  theory,  many  new  random  variables  were 
included  and  various  distribution  functions  studied  in  order 
to  understand  better  this  type  of  casualty  occurence 
paradigm. 
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This  thesis  only  dealt  with  the  case  of  homogeneous 
forces  which  were  using  the  same  acquisition  mode  to  attrit 
tne  targets,  which  were  passive  X  forces.  Thus  we  could  use 
an  ordinary  renewal  process  to  study  the  behaviour  of  casu- 
alty occurencas.  We  emphasized  the  probabilistic  function 
of  each  r.v.,  so  once  we  found  the  distribution  of  interar- 
rival  times  of  casualties,  casualty  analysis  could  be  done 
by  plugging  distribution  functions  into  the  eguation  devel- 
oped. Finding  the  distribution  of  interarrival  times  of 
casualties  was  the  critical  factor  in  continuing  the  casu- 
alty   analysis. 

The  final  conclusion  is  that  total  casualty  occurence  in 
co§bat  can  be  assumed  to  be  a  Poisson  process ,  and  thus 
tiroes  between  casualties  are  exponentially  distributed  if  we 
assume   a  large    number  of  combatants. 

We  have  dealt  here  only  with  a  passive  X  f orce  (  i.e.  one 
that  does  not  shoot  back  )  ,  perhaps  an  unrealistic  assump- 
tion. But  still  this  study  permits  tetter  understanding  of 
combat  attrition,  through  the  use  of  probability  theory, 
than    was    previously    available. 

B.       SUGGESTIONS    FOR    FDETHEE    STUDY 

In  this  thesis,  we  have  only  dealt  with  passive 
targets.  A  more  complex  formulation  may  be  required  for  real 
combat,  since  enemy  targets  shoot  back  and  cause  firer 
attrition. 

In  addition,  a  delayed  renewal  process  could  be  applied 
for  the  case  where  different  acquisition  modes  are  used  for 
the  first  shooting  and  for  each  subsequent  shootings.  A 
reinforcement  of  the  forces  could  also  be  included. 
Casualty  analysis  could  then  be  done  using  the  theory  of 
alternating  renewal  processes,  once  we  know  the  rate  of 
reinforcement. 
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APPENDIX    A 
FORTRAN    PROGRAMMING   FOR    SIMULATING    THE    KILL   TIME     (N    =    3) 


REAL    01 (3)  ,02(3] £U3, U 4 - T3IN, TIME  13)  .INTIME (40, 14) 
INTEGER    JIGI,FINDEX,L,H#  N,KIGT,  IGT  (14)  ,RTGT  (14) 
DOUBLE   PRECISION    DSEED 
ESEED    =    123457. CO    C 
DO    200   K    =    1,40 

M    =    14 

N    =    3 

I    a     1     C 

CALL   3GUSS (DSEED, N, 01)     C 
DO     100    I    =     1,N 
TIME  (I)     =     01(1) 
100  CONTINUE    C 

DO     105    J    =     1  ,M 
RTd  (J)    =    J 
105  CONTINUE    C 

DSEED   =     DSEED    +    23.  DO    C 
CALL   3GUBS (DSEED, N,  02)     C 
DO     110    I    =     1,N 

TGT(I)     =    INT(U2(I)     *    FLOAT(H))     +     1 
110  CONTINUE    C 

10  CONTINUE 

TMIN   =    2. 
CO    130    I    =    1,N 

IF(IMIN.LE.1IME(I))     GO    TO    130 
TMIN    =    TIME  (I) 
KTST    =    TGT  (I) 
FIND EX    =    I 
130  CONTINOE    C 

INTIME (K,L)      =    TMIN    C 
DO    125    J    =    1  ,M 

I?(STGT(J)  .NE.KTGT)     GO    TO    125 
JTGI    =    J 

RTGI  (JTGT)     =    RTGT  (M) 
125  CONTINUE    C 

H    =    M    -     1 
L    =    L    +     1    C 
DO    135    I    =    1,N 

IF(KTGT.EQ.IGT  (I)  )     GO    TO    70 
TIME (I)     =     TIME    I       -    IMIN 
GC   TO     135 
70       '  DSEED    =    DSEED    +    1 2. DO 

CALL    GGU3S  (DSEED,  1,03) 
TIME (I)    =    03 
DSZED    =    DSEED    +    3  .  DO 
CALL     GGU3S  (DSEED,  1,04) 
TGT(I)     =    HIGT{INT  (U4*FL0AT  (M)  )  +1) 
135  CONTINUE    C 

IFfM.GE.  1)     GC    TO    10 
20  0       CONTINUE 

WRITE  (6,300)  (  (INTIME  (X,L)  ,  L=  1  ,  14)  ,K=1,40) 
STOP 
300       F0£MAT(1X,  14F5.2) 
END 
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APPENDIX    B 
SIMULATED   DATA   FOE    KILL    TIME    (N    =    3,    M   =    6) 


0. 

260  7 

0. 

5056 

0.2000 

0. 

4392 

0. 

1057 

0. 

.0649 

0. 

3914 

0. 

3860 

0.1659 

0. 

3315 

0. 

0303 

0. 

.1646 

0. 

1156 

0. 

2987 

0.01 87 

0. 

3911 

0. 

0438 

0. 

,4744 

0. 

4545 

0. 

0154 

0.13  12 

0. 

2905 

0. 

1360 

0. 

,3572 

0. 

4138 

0. 

1269 

0.0678 

0. 

,1355 

0. 

1724 

0. 

.0101 

0. 

088  4 

0. 

0730 

0.2945 

0. 

0234 

0. 

7305 

0. 

.08  20 

0. 

4117 

0. 

3926 

0.0363 

0. 

0o97 

0. 

0277 

0. 

.0151 

0. 

0499 

Q. 

2410 

0.3564 

0. 

0033 

0. 

1442 

0. 

.2979 

0. 

4637 

0. 

2188 

0.1503 

0. 

,054o 

0. 

5555 

0. 

.0280 

0. 

7428 

c. 

0736 

0.3629 

0. 

0375 

0. 

,0319 

0. 

.0322 

0. 

232  6 

0. 

3066 

0.2037 

0. 

2932 

0. 

1050 

0. 

.1466 

0. 

1105 

0. 

5015 

0.2316 

0. 

1160 

0. 

,2916 

0. 

.4032 

0. 

0943 

0. 

1942 

0.2346 

0. 

0342 

0. 

2771 

0. 

.0218 

0. 

5372 

0. 

0104 

0.3656 

0. 

143  1 

0. 

0724 

0. 

.1435 

0. 

0492 

0. 

0554 

0.1336 

0. 

1750 

0. 

5712 

0. 

.0609 

0. 

2496 

0. 

0227 

0.11 02 

0. 

2093 

3. 

2394 

0. 

.1294 

0. 

453  0 

0. 

0027 

0.1524 

0. 

0436 

0. 

0132 

0. 

.6276 

0. 

2916 

0. 

2342 

0.0265 

0. 

,5775 

0. 

,3016 

0. 

.1875 

0. 

229  9 

0. 

0756 

0.5040 

0. 

1737 

0. 

0637 

0. 

.5797 

0. 

7527 

0. 

0343 

0.1628 

0. 

0063 

0. 

3606 

0. 

.0377 

0. 

220  1 

0. 

3  036 

0.1323 

0. 

5570 

0, 

,0256 

0. 

.3332 

0. 

6097 

0. 

0538 

0.05  16 

0. 

3331 

0. 

1933 

0, 

.0139 

0. 

782  3 

0. 

0719 

0.1470 

0. 

,01ob 

0. 

0743 

0. 

.0118 

0. 

1133 

0. 

6589 

0.0144 

0. 

0712 

0. 

,  1216 

0, 

.2234 

0. 

003  5 

0. 

0602 

0.44  10 

0. 

2291 

0. 

1626 

0. 

.2609 

0. 

0174 

0. 

56  69 

0.1626 

0. 

2062 

0. 

,  1578 

0, 

.0334 

0. 

6032 

0. 

00C3 

0.25  92 

0. 

0024 

0. 

0723 

0. 

.2653 

0. 

3186 

0. 

04  92 

0.39  18 

0. 

1603 

0. 

,  1857 

0, 

.0369 

0. 

556  4 

0. 

0105 

0.0674 

0. 

,2037 

0. 

0184 

0. 

.1233 

0. 

570  6 

0. 

3844 

0.3504 

0. 

2022 

0. 

0417 

0. 

.1048 

0. 

2102 

0. 

4294 

0.11 55 

0. 

0393 

0. 

,  1705 

0, 

.4165 

0. 

019  1 

0. 

2929 

0.4283 

0. 

2152 

0, 

,  1427 

0, 

.4954 

0. 

1570 

0. 

0601 

0.31  1  1 

0. 

4439 

0. 

0499 

0. 

.0377 

0. 

193  1 

0. 

4658 

0.0251 

0. 

15b  8 

0. 

,0375 

0 

.2395 

0. 

1372 

0. 

41  83 

0.1055 

0. 

2174 

0. 

,0335 

0. 

.2340 

0. 

0523 

0. 

2950 

0.4356 

;>. 

,0230 

0. 

1054 

0, 

.0926 

0. 

302  4 

0. 

5260 

0.0669 

0. 

,0079 

0. 

,  1632 

0. 

.0743 

0. 

0645 

0. 

1  504 

0.0088 

0. 

0741 

0. 

5581 

0, 

.0223 

0. 

0220 

0. 

4350 

0.2320 

0. 

,0207 

0. 

,2634 

0, 

.2070 

0. 

126  5 

0. 

1379 

0.1063 

0. 

.0695 

0, 

.0050 

0, 

.5903 
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APPENDIX    C 
DISCRETE    APPROXIMATION    FO    RENEWAL    FUNCTION 


SJCB 

1  REAL    P0,M1  (401)  ,K  (41)  ,IIME  (41) 

2  M1 (1)     =    0. 

3  PO    =     1./100- 

4  DO    100   I    =  1.400 

5  IF  (I.GT.1  00)  GO    TO    200 

6  H1  (1+1)    »  H1  (I)     +     (1.     +    M1(I))     *    PO 

7  GO    IO    100 

8  200            M  1(1+1)     ■  M1  (I)     ♦     (1.     +    M1(I))     *    PO 

*      -     h-    *    M1  (1-100)  )  *    PO 

9  100       CONTINUE 

10  DO     1  10    J    =  1.4  1 

11  TIME(J)     =  0.  1    *     (J-1) 

12                      a  (j)   =  in  (1  +  10*  (j-1) 

13  WRITE  (6,*)     TIME(J)  ,21  (J) 

14  110       CONTINUE 

15  STOP 

16  END  0  SENTRY 
0.0000000  0.0000000 
0.  1000000  0.  1046219 
0.2000000  0.2201895 
0.3000001  0.3478478 
0.4000001  0.4883622 
0.5000001  0.6446293 
0.6000001  0.3166940 
0.7000002  1.0067590 
0.8000002  1.2167040 
0.9000002  1.4486160 
1.0000000  1.7047900 
1.0999990  1.8733930 
1.1999990  2.0587170 
1.3000000  2.2452o50 
1.3999990  2.4373700 
1.5000000  2.6341510 
1.6000000  2.3344850 
1.6999990  3.0369570 
1.3000000  3.2393290 
1.8999990  3.4409610 
2.0000000  3.6377740 
2.  1000000  3.  3320340 

2.  1999990  4.0281400 
2.3000000  4.2256030 
2.4000000  4.4239440 
2.5000000  4.6227140 
2.Q000000  4.8215160 
2.6999990  5.0200490 
2.8000000  5.2181320 
2.9000000  5.4157760 
3.0000000  5.6132290 

3.  1000000  5.8109340 
3.2000000  o. 0089100 
3.3000000  6.2070170 
3.4000000  6.4051470 
3.5000000  6.6032300 
3.6000000  6.8012370 
3.7000000  6.9991710 
3.3000000  7.1970o00 
3.9000000  7.3949500 
4.0000000  7.5928790 
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